OPEN 3 ACCESS Freely available online 



<P» PLOS I ONE 



A Method for Quantitative Analysis of Standard and (S\ 
High-Throughput qPCR Expression Data Based on Input &SL 
Sample Quantity 

Mateusz G. Adamski 1,2 , Patryk Gumann 3,4 , Alison E. Baird 1 * 

1 Neurology, SUNY Downstate Medical Center, Brooklyn, New York, United States of America, 2 Neurology, UJCM, Krakow, Poland, 3 Department of Physics, Harvard 
University, Cambridge, Massachusetts, United States of America, 4 Institute for Quantum Computing, Department of Physics and Astronomy, University of Waterloo, 
Waterloo, Ontario, Canada 



Abstract 

Over the past decade rapid advances have occurred in the understanding of RNA expression and its regulation. Quantitative 
polymerase chain reactions (qPCR) have become the gold standard for quantifying gene expression. Microfluidic next 
generation, high throughput qPCR now permits the detection of transcript copy number in thousands of reactions 
simultaneously, dramatically increasing the sensitivity over standard qPCR. Here we present a gene expression analysis 
method applicable to both standard polymerase chain reactions (qPCR) and high throughput qPCR. This technique is 
adjusted to the input sample quantity (e.g., the number of cells) and is independent of control gene expression. It is 
efficiency-corrected and with the use of a universal reference sample (commercial complementary DNA (cDNA)) permits the 
normalization of results between different batches and between different instruments - regardless of potential differences 
in transcript amplification efficiency. Modifications of the input quantity method include (1) the achievement of absolute 
quantification and (2) a non-efficiency corrected analysis. When compared to other commonly used algorithms the input 
quantity method proved to be valid. This method is of particular value for clinical studies of whole blood and circulating 
leukocytes where cell counts are readily available. 
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Introduction 

Over the past decade a rapid increase has occurred in the 
understanding of RNA expression and its regulation. Quantitative 
polymerase chain reaction(s) (qPCR) have become the gold 
standard for measuring gene expression. Accurate analysis of 
qPCR data is crucial for optimal results and a number of well- 
defined methods are in use to calculate gene expression. These 
include the comparative C x method [1], the efficiency corrected 
method [2] and sigmoidal curve fitting methods [3] , all of which 
provide relative quantitative information. A standard curve of 
serial dilutions of a known sample is additionally required to 
measure the absolute number of transcript copies in a sample. 

For most scientific purposes, relative quantification, expressed as 
fold change, is sufficient to provide the required information. 
Hence, the comparative Ct and efficiency corrected methods, as 
well as the sigmoidal curve fitting methods are widely employed, 
but each method has strengths and weaknesses. The comparative 
C T method by Livak et al. [1] has the advantage of ease of use but 
is based on the assumption that transcript amplification efficiencies 
are 100%. In the efficiency corrected method by Pfaffl [2] the 
relative expression ratio is calculated only from the real-time PCR 
efficiencies and the crossing point deviation of an unknown sample 
versus a control. This model needs no calibration curve and gives 



improved quantification but is complex to use and requires 
determination of the amplification efficiency. 

Furthermore, all of these methods require the use of reference 
(control or housekeeping) genes to correct for unequal amounts of 
biological material that may exist between the tested samples. The 
commonly used housekeeping genes were initially selected on the 
basis of their abundance and expression in a wide variety of tissues. 
An absolute requirement and widely held assumption of house- 
keeping genes has been that their expression is constant under all 
conditions and is unaffected by the experimental conditions [4]. 
However, the expression of commonly used housekeeping genes 
has since been found to vary considerably in many conditions [5— 
12]. In the case of in vitro or ex-vivo experiments it is usually 
possible to perform additional experiments to identify and validate 
appropriate control genes. In the case of clinical studies, however, 
where sample volumes are usually limited, it is rarely possible to 
test gene expression before and after the experiment (i.e., before 
and after the disease occurs). 

The advent of next generation high throughput qPCR, based 
on reaction volumes scaled to the nanoliter range and with a 
consequent dramatic reduction in the volume of reagents and 
samples, has been a major advance for the analysis of clinical 
samples [13]. The Fluidigm Biomark system, one of the new high- 
throughput reverse transcription PCR (HT RT-qPCR) systems, 
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permits up to 96 transcripts in 96 samples to be studied 
simultaneously during a single run, in a total of 9216 reactions. 
This allows many more transcripts to be studied from routine 
clinical samples, representing a 40 to 50 fold improvement in 
efficiency over standard qPCR [14,15]. However, HT RT-qPCR 
has also raised new issues; for example, transcript amplification 
efficiency may be affected by potential interactions (i.e., primer 
dimer, competition) between multiple primers during the pre- 
amplification and amplification steps. 

Here, we present a method for the measurement of the absolute 
gene expression for standard and high throughput qPCR 
experiments based on the input sample quantity. Based on this 
method three equations were developed: (1) for the measurement 
of fold change differences between target and control samples; (2) 
for the comparison of results from different experiments and 
different machines after normalization to a reference cDNA 
sample; (3) for analyses of samples of unknown efficiency. Gene 
expression results calculated using the input quantity method were 
then validated in a serial dilution series of commercial cDNA and 
using different starting cell concentrations. In clinical samples, fold 
change values calculated with the input quantity method were 
compared to values obtained using other commonly used 
algorithms. The input quantity method has the advantages of 
avoiding the use of control genes, of being efficiency corrected, 
and providing both fold change and absolute results. This method 
can also be applied in the verification and quantification of 
qualitative results from microarray studies for multiple genes. 

Theory 

1. Requirements for the input quantity method 

The input quantity method has several requirements. First, the 
amount of material used for RNA extraction has to be measured: 
for example, cell count is required for cell suspensions (e.g., 
peripheral blood mononuclear cells (PBMCs), lymphocytes and 
cell lines), white blood cell (WBC) counts are needed for whole 
blood studies and tissue volumes are needed for solid tissues. 
Secondly, for reverse transcription of RNA to cDNA the same 
reagents, volumes and protocols for a given experiment need to be 
used. Thirdly, the amplification efficiency and correlation 
coefficients (R 2 ) should be assessed for each gene assay based on 
a standard dilution series. Finally, full application of this method 
requires the use of a standard sample (i.e., commercial cDNA - 
reverse transcribed cDNA from RNA extracted from all human 
tissues) for each measurement. 

2. Mathematical model for qPCR amplification 

As per Livak et al. [1], in the qPCR target cDNA sequence is 
amplified in an exponential fashion: 

J„=J„x(l + £)" (1) 

where X n is the number of target cDNA molecules after n cycles, 
X 0 is the number of cDNA molecules before amplification, E is the 
efficiency of target cDNA amplification and n is the number of 
amplification cycles. In the case of perfect efficiency (E = 100%) 
the number of target cDNA molecules doubles every cycle. 

In qPCR, the number of target cDNA molecules for a given 
sample is reflected by the threshold cycle - or according to the 
MIQE guidelines [4], quantification cycle (Cq) - because Cq is the 
intersection between an amplification curve and threshold. The 
threshold is the level of fluorescence above background fluores- 
cence - set at the same level for all samples in the experiment. 
Each sample that crosses the threshold (regardless of the 



amplification cycle number) has the same fluorescence intensity 
hence the same target cDNA copy number. 

X nCq = X 0 x(l+E)" c « = K (2) 

where XnCq is the number of target cDNA molecules at the Cq, 
nCq is the cycle number at which amplification crosses the 
threshold and K is a constant value for all samples in a given 
experiment. 

3. Analysis normalized to input sample quantity 

In order to adjust the results of gene expression to unequal 
amounts of starting material the number of cells used for RNA 
extraction has to be incorporated into Equation 2. 

X 0 = X c xcc (3) 

where Xc is the transcript number per cell and cc is the number of 
cells used for RNA extraction (e.g., complete blood count for 
whole blood analysis, or hemocytometer cell count for cell subset 
analysis). Hence, 



K = (X c xcc)x(l + E) nCq (4) 

Therefore to compare gene expression between target (T) and 
control (C) samples where E and K are the same for T and C, ccT 
is the input cell count for target sample and ccC is the cell input for 
the control sample. For the target samples the following formula is 
obtained: 

K = (T c xccT)x{\+E) nCcl ' T (5) 

where Tc is the number of transcripts per cell in the target 
samples. 

For the reference or control samples the following formula is 
obtained. 

K = (C c xccC)x{\+E) nCq - c (6) 

where Cc is the number of transcripts per cell in the reference 
samples. 

As K is constant, Equations 4 and 5 equal each other: 

(T c x ccT) x (1 +Ef c i' T = (C c x ccC) x (1 +E)" Cc '' c (7) 

To obtain the comparison between target and control samples: 
g = g x( l +jE) («C ? ,C-«C ? ,71 (8) 

This way we can obtain the measure of gene expression 
expressed as a fold change difference between the test and control 
samples. 

4. Analysis normalized to input quantity and normalized 
to standard cDNA 

When a standard reference sample is introduced, for example a 
sample that contains a high concentration of studied transcripts, 
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the following modifications are made, starting with Equation 2, K 
for sample X with a starting quantity of cc is: 



K = (X c xcc)x(l+E) nC ' l - x 



(9) 



K for a standard cDNA of uniform quantity is: 

K = cDNA 0 x (1 + E) nCq ' cDNA (10) 

Normalizing to cDNA: 
(X c x ccX) x (1 + E)" Cq ' x = cDNAo x (1 + Ef c "' cDNA (11) 



X c = cDNA 0 x 1 - ; ^ (12) 



Since the number of transcripts before amplification in standard 
cDNA (cDNAq) is constant we may assume it is equal to 1 then: 



Xc. 



(l+E)< 



(nCq,cDNA-nCq,X) 



CCX 



(13) 



To obtain the comparison between test and control samples, the 
respective T r and C r are calculated using Equation 13. Then T c is 
divided by C c to obtain the measure of gene expression, expressed 
as a fold change. 

5. Analysis normalized to input quantity and/or 
normalized to standard cDNA without known efficiency 

If E for the working primers is not assessed in the experiment, 
one may make an assumption that the E equals 100% -then 
Equation 8 is: 



7o _ cc,C 
C 0 cc,T 



. 2(nCq,C-nCq,T) 



(14) 



Whereas, adjusting to the standard cDNA sample, for sample X 
Equation 12 is: 



2(nCq,cDNA-nCq,X) 

ccC 



(15) 



Materials and Methods 

To assess the reliability of the input quantity method, the 
stability of expression values calculated across serial dilutions of a 
standard cDNA sample and of different starting numbers of two 
samples of peripheral blood mononuclear cells (PBMCs) were 
determined. The validity of the input quantity method was 
assessed by comparison to fold changes obtained using the Livak 
[1] and Pfaffl [2] methods for three transcripts in a cohort of stroke 
patients and control subjects. 

The Institutional Review Board at the State University of New 
York (SUNY) Downstate Medical Center approved the study. All 



study participants and/ or authorized representatives gave full and 
signed informed consent. Where applicable, the conduct and 
reporting of the study are in accordance with the MIQE criteria 
[4]. The detailed laboratory protocols but not the data analysis 
described in this manuscript have been previously published [14]. 

1. RNA extraction and reverse transcription 

Whole blood was obtained from 38 ischemic stroke patients 
between 7 and 90 days post stroke and from 17 sex- and race- 
matched control subjects. RNA was extracted using column 
separation (All-in-One Kit; Norgen Biotek, Thorold, Ontario, 
Canada) from 100 |Ol of whole blood and from a median of 2.0 
million CD4 + cells. Peripheral blood mononuclear cells (PBMCs) 
from two control subjects were used for the cell dilution 
experiment, with RNA isolated from triplicate samples of 2 
million, 1 million, 0.5 million and 0.25 million cells. Cellular 
counts (millions of cells per |ll) were measured using a 
hemocytometer for CD4 + and for PBMCs; for whole blood, the 
total white blood cell count was obtained from the laboratory- 
measured complete blood count (CBC) in each study subject. 

Density gradient centrifugation with Histopaque 1077 and 1119 
(Sigma-Aldrich, St. Louis, MO) was used to separate the PBMC 
fraction from the whole blood. Positive magnetic bead separation 
(Miltenyi Biotec, Bergisch Gladbach, Germany) was used to 
separate CD4 + from PBMCs - the cellular purity was over 97%. 
The extracted RNA was resuspended in 50 JLll of elution solution 
(All-in-One Kit protocol). cDNA was synthetized using the High 
Capacity cDNA Reverse Transcription Kit (Life Technologies, 
Carlsbad, CA), based on random hexamers, according to the 
manufacturer's protocol. Following the protocol, the proportion of 
RNA solution to 2x RT master mix was 1:1. 

2. Primer development, RT qPCR and HT-RT qPCR 

The primers for qPCR were self-designed, commercially 
synthesized by Invitrogen and wet tested using standard RT 
qPCR (StepOnePlus Real-Time PCR Systems; Applied Biosys- 
tems). 

Standard RT qPCR (StepOnePlus Real-Time PCR Systems; 
Applied Biosystems) was used to measure the expression of FDFT1 
in the cell dilution experiment. Each sample and no template 
control were measured in triplicate. Based on a standard dilution 
series the efficiency for FDFT1 in this experiment was 94%. 

HT RT-qPCR was run on the BioMark HD System, using 
96x96 Fluidigm Dynamic Arrays (Fluidigm, South San Francisco, 
CA). HT-RT qPCR was used first, to measure the expression of 
FUT4, CD3E, FDFT1 and B2M in serial dilutions of commercial 
cDNA (Universal cDNA Reverse Transcribed by Random 
Hexamer: Human Normal Tissues; Biochain, Newark, CA) and 
second, to compare the expression of FDFT1, CD3E and B2M 
between control subjects and stroke patients in whole blood and 
CD4 + T lymphocytes. Two 5 point, four-fold serial dilution series 
of commercial cDNA were run in triplicate on two different plates. 
The volumes of commercial cDNA (diluent) in each dilution were: 
100 m_1 (1:1), 25 m_1 (1:4), 6.25 |ll (1:16), 1.5625 ul (1:64) and 
0.39 |ll (1:256). According to the manufacturer's protocol, the 
assay for each HT RT-qPCR experiment contained 10 |0l of 
cDNA. The efficiencies for the genes, assessed with HT RT- 
qPCR, were: B2M- 87%, FDFT1- 86%, FUT4- 79% and CD3E- 
79%. Five separate gene expression plates were used in this 
experiment. To normalize the gene expression results for stroke 
and control samples from different plates, a sample of commercial 
cDNA (containing high concentrations of all of the transcripts 
studied) of standard concentration and volume was run in 
duplicate on each plate. Each raw gene expression result 
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(expressed as Cq) was normalized to the average Cq value for the 
same gene in the commercial cDNA samples that were run on the 
same plate (sample Cq value for gene X was subtracted from the 
average commercial cDNA Cq for gene X). 

3. Calculation of fold changes 

Fold change differences between stroke patients and control 
subjects for B2M and CD3E were calculated using the input 
sample quantity method according to Equation 13. The relative 
gene expression for B2M and CD3E were measured using the 
comparative Cx method of Livak et al. [1] and the efficiency 
corrected method of Pfaffl [2] . For these calculations FDFT1 was 
used as control gene as its expression was not different in stroke 
patients compared to control subjects, based on the input quantity 
method (p>0.05). 

4. Statistical analyses 

The statistical analyses were performed using "R", version 
2.15.2. For the cDNA dilution analysis, linear regression modeling 
was used. For the cell dilution series, the data were analyzed using 
one way ANOVA, Welch's correction for inhomogeneity of 
variances and post hoc t. tests with false discovery rate correction. 
For the analysis of the stroke versus control data, the 95% CI for 
the fold change values were calculated using the R package 
"mratios" and Dunnetts method; Wilcoxon rank sum tests were 
used for between group comparisons. 

Results 

1. Gene expression measurements across different input 
volumes of a standard cDNA sample 

To confirm the reliability of the sample input quantity method 
the expression of 4 transcripts (FUT4, CD3E, FDFT1 and B2M) 
was measured in 5 point and 4-fold two serial dilutions of a 
standard cDNA sample. To measure the concentrations of each of 
the four transcripts in the standard cDNA sample, the results were 
normalized to the volume of diluent 100 Ml (1:1), 25 ul (1:4), 
6.25 ul (1:16), 1.5625 pj (1:64) and 0.39 u.1 (1:256). Using this 
normalization procedure the same expression values were 
expected across the range of dilutions of the standard cDNA 
sample. The samples were run in triplicate on two separate plates 
giving 6 readings per input volume. The expression of all four 
genes calculated with the input quantity method was stable 
(Table 1, Figure 1). Detailed data are provided in a supplemental 
table (Table SI). 

2. Reliability of gene expression measurements across 
different starting numbers of cells 

In order to determine the influence of variables present prior to 
the RT qPCR step (cell counting, RNA isolation and RT PCR) 
the expression of FDFF1 in different starting numbers of PBMCs 
from two control subjects was measured. The raw data were 



normalized to the starting number of cells for each subject. The 
starting numbers of cells (2 million, 1 million, 0.5 million and 0.25 
million) were within the range of the manufacturer's recommen- 
dations for RNA extraction (All-in-One Kit, Norgen Biotec). 

Based on the input quantity method the expression of FDFT1 
was significandy different across the input cell counts for both 
subjects (p = 1.4e-7, Subject 1 and p = 5.5e-5, Subject 2) (Table 2). 
Post hoc tests revealed that the expression of FDFT1 in the 0.25 
million input cell count in both subjects differed significantly from 
the other input cell concentrations: in Subject 1 (versus 2 million, 
p = 2.7e-6, versus 1 million, p = 0.000 16 and versus 0.5 million, 
p = 7.6e-5) and in Subject 2 (versus 2 million, p = 5.9e-5, versus 1 
million, p= 1.3e-6 and versus 0.5 million, p= 1.7e-6). Compari- 
sons between the 2 million, 1 million and 0.5 million input cell 
numbers were not statistically significant for both subjects (p< 
0.05). Detailed data are provided in a supplemental table (Table 
S2). 

3. Expression of CD3E and B2M in the late phase of stroke 
and in control subjects calculated using three methods 

To assess the validity of the input quantity method using clinical 
samples, the expression of CD3E and B2M in whole blood and in 
CD4 + T lymphocytes was compared between patients in the 
delayed phase of stroke and control subjects. Fold change 
differences in gene expression were measured using the input 
quantity method (normalized to cell count), and the Livak and 
Pfaffl methods. 

By all methods B2M expression was significantly increased in 
whole blood in the delayed phase of stroke and CD3E was 
significantly increased in CD4 cells (Table 3). No alterations in the 
expression of CD3E were found in whole blood. A borderline 
increased in B2M expression in CD4 cells was found using the 
input quantity method. Detailed data are provided in a 
supplemental table (Table S3). 

Discussion 

Several gene expression analysis methods are in common use, 
but the input quantity approach presented here offers two major 
advantages. Firstly, this method is independent of control genes. 
Secondly, with the assumptions of 1) uniform efficiency of RNA 
extraction and RT qPCR and 2) a constant concentration and 
volume of a standard sample, this method permits absolute 
quantification, expressed as the fraction of transcripts in the 
standard sample, across different experiments. The proposed 
algorithm is efficiency corrected, although analysis of results 
without known efficiency is also possible. With the use of a 
standard sample, the input quantity method also permits the 
comparison and analysis of results from different batches and 
results acquired on different qPCR machines. Furthermore, with 
the advent of HT RT-qPCR, this analytical method is also very 
useful for clinical research, where sample volumes are limited. 



Table 1. Expression of FUT4, CD3E, FDFT1 and B2M across serial volumes of a standard cDNA sample. 







FUT4 


CD3E 


FDFT1 


B2M 


Coefficient 


-4.3e-8 


-6.43e-8 


-3.8e-7 


-3.3e-6 


P value 


0.49 


0.65 


0.61 


0.48 


R 2 


-0.018 


-0.028 


-0.026 


-0.016 



Dilution coefficient, p and R 2 values were obtained from linear regression analysis for each transcript. 
doi:1 0.1 371 /joumal.pone.01 0391 7.t001 
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Figure 1 . Expression of FUT4, CD3E, FDFT1 and B2M'm a standard dilution series of reference cDNA sample normalized to the volume 
of diluent using sample input quantity method. 

doi:1 0.1 371 /journal.pone.01 0391 7.g001 



Table 2. Expression of FDFT1 in cell dilution series. 







2 million cells 


1 million cells 


0.5 million cells 


0.25 million cells 


P 


Subject 1 


0.26±0.01 


0.23 ±0.02 


0.24±0.06 


0.1 5 ±0.02** 


«0.01 


Subject 2 


0.049±0.003 


0.041 ±0.005 


0.043±0.015 


0.072±0.013** 


«0.01 



p values were calculated using a one-way ANOVA. 

**Post hoc tests revealed that expression of FDFT1 in the 0.25 million input cell count differed significantly from the other input cell concentrations in both subjects. 
doi:1 0.1 371 /journal.pone.01 0391 7.t002 
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Table 3. Fold change difference in the expression of B2M and CD3E in late phase stroke versus control subjects. 





B2M 


B2M 


B2M 


CD3E 


CD3E 


CD3E 


Whole 
blood 


Input Quantity 
Method 


Livak 


Pfaffl 


Input Quantity 
Method 


Livak 


Pfaffl 


Fold Change 


2.51 


2.19 


2.28 


1.27 


1.12 


1.22 


95% CI 


1.26, 15.89 


1.26, 5.94 


1.32, 6.20 


0.67, 3.34 


0.70, 2.01 


0.79, 2.07 


P 


0.017 


0.006 


0.003 


0.19 


0.48 


0.42 


CD4 


Fold Change 


1.35 


0.57 


0.70 


3.13 


1.78 


2.10 


95% CI 


0.94, 2.17 


0.26, 1.15 


0.41, 1.23 


1.61, 25.8 


1.16, 3.42 


1.35, 4.25 


P 


0.02 


0.4999 


0.26 


2.10e-05 


0.0084 


7.50e-05 



Fold change and 95% CI calculated using the 3 methods; data analyzed using R package "mratios" and Dunnetts method; Wilcoxon rank sum tests used for between 

group comparisons. 

doi:1 0.1 371 /journal.pone.01 0391 7.t003 



Our analyses show that the sample input quantity method 
permits gene expression to be measured across a wide range of 
commercial cDNA. Although the performance of both RNA 
extraction and RT qPCR may differ significantly across different 
cell concentrations and kits [15], our results show that, using the 
same protocol and reagents within the input quantities we tested, 
these variables can be successfully controlled. Furthermore, the 
expression of B2M and CD3E in study subjects calculated using 
three methods was highly concordant. 

The rationale for the use of housekeeping (or control or 
reference genes) is to correct gene expression results, reflected as 
differences in Cq values between target and control samples, that 
could result from two main factors: different amounts of starting 
material or different levels of expression. Traditionally, house- 
keeping genes have been chosen on the basis of their abundance, 
ubiquitous expression across tissues and the assumption that their 
expression is stable under physiological and experimental condi- 
tions. However, the expression of conventionally used housekeep- 
ing genes varies considerably in many conditions. Therefore, 
reference gene selection requires additional experiments to 
validate gene expression stability under different experimental 
conditions [6-12,14]. In many conditions, especially in the clinical 
setting, it is not possible to measure the effect of the disease/ 
condition on reference gene expression. 

The algorithm used for our sample input quantity method 
employs normalization to the sample input quantity (cell count, 
tissue volume etc.), which in result permits an absolute gene 
expression analysis. This method varies from the relative analysis 
approach, where results are normalized to reference gene 
expression. Due to normalization to the input quantity (measured 
in absolute scale) the measure of gene expression remains absolute, 
as in our method. In contrast, the gene expression from the 
relative analysis approach is based on the normalization to 
reference gene expression. Thus the ratio of the target gene 
expression to the reference gene expression represents a relative 
measure. By introducing a standard sample (of a stable transcript 
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concentration), our method allows us to compare gene expression 
between different experiments. Instead of directly measuring 
transcript copy number- as it is commonly done in absolute 
measurements of gene expression- in our method, the measured 
gene expression is presented as a fraction of transcripts present in 
the standard sample. This fraction can be converted to the 
transcript copy number by measuring concentration of the target 
gene in the standard sample. 

The input quantity approach presented here can be applied to 
clinical studies, to verify and quantitate microarray results, and to 
large scale studies of gene or microRNA expression. Having 
knowledge of the input cell count for all samples and the use of a 
uniform standard, first, allows normalization to the amount of 
starting material, and second, the use of the same standard allows 
normalization of results between different laboratories and 
different equipment. 
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